HIP - 1997 - 49 / TH 



Low-Energy Pion-Nucleon Interaction and the Sigma- Term* 

M.E. Sainio 

Department of Physics, University of Helsinki, P.O. Box 9, OOOI4 Helsinki, Finland 

Abstract 

A dispersion framework, appropriate for discussing the low-energy pion-nucleon in- 
teraction, is reviewed. Sensitivity of the isoscalar ZJ-ampUtude at the unphysical Cheng- 
Dashen point to input from different phase shift analyses and low-energy experiments is 

- - - , discussed. 
!>■ 

Q> . INTRODUCTION 

In the limit of vanishing u- and d-quark masses the QCD Hamiltonian has chiral SU{2) x 
SU (2) symmetry. The symmetry group would be SU{'i) x SU{2>), if also the strange quark is 
• treated as massless. The QCD vacuum does not, however, have the same symmetry, and the 

symmetry is said to be hidden, or spontaneously broken. As a consequence, the Goldstone 
theorem dictates that massless pseudoscalar mesons appear. They acquire a mass through 
the small, but nonzero, quark masses. The quark masses also shift the mass of the proton. 
The proton matrix element of the light quark mass term in the Hamiltonian, the sigma-term, 
, is a measure of the explicit chiral symmetry breaking. 
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0^ ; a = —{p\uu + dd\p) , (1) 
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where rh = |(m„ -|- md), m is the proton mass and \p) denotes the physical one-proton state 
Qh! normalized as {p'\p) = (27r)^ 2po 6{p' — p). Algebraically a can be written in the form 
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THE SIGMA TERM 

The sigma-term is defined as 



TTi {p\uu + dd — 2ss\p) . 
^~2^ l-y ' ^' 

where the parameter y, the strange quark content of the proton, is defined as 

{p\uu + dd\p) 

The combination cj(l — y) can be calculated in Chiral Perturbation Theory (CHPT), and 
with ms/m ~ 25 an estimate can be given in different orders of the quark mass expansion 

o-(l-y) ~ 26 MeV 0{mg) 
~ 35 ±5 MeV C'("^g^^) 

~ 36 ±7 MeV 0{ml). (4) 

The corrections to the leading order calculation refiect the meson cloud effects. The result of 

the order mq^"^ is from Ref. [1,2] and the result from Ref. [3]. In the OZI rule limit {y = 0) 
a value for a follows to a given order of the quark mass expansion. 

* Invited talk presented at the "7th International Symposium on Meson-Nucleon Physics and the Structure of 
the Nucleon", Vancouver, B.C., Canada, July 28- Aug. 1, 1997. 
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A low-energy theorem of chiral symmetry allows for another determination of the cj-term. 
Namely, the pion-nucleon scattering amplitude at a particular unphysical (but on-shell) point, 
so called the Cheng-Dashen point, can be related to the nucleon scalar form factor up to 
corrections of order log ji^ [4] 

S = F^D+{u = 0,t = 2f/) = a{2fi^) + Ar. (5) 

Here = 92.4 MeV is the pion decay constant, Z)+ denotes the isoscalar combination of the 
pion-nucleon amplitudes D = A+uB, where the pscudovcctor Born term has been subtracted, 
i.e. D = D — Dp^, u = {s — u)/Am in terms of the Mandelstam variables and /x is the pion 
mass. The scalar form factor a{t) is defined by 

u{p')a{t)u{p) = {p'\m{uu + dd)\p) ; t = {p'—p)'^. (6) 

The correction term, A/j, is formally of order fi^, and to one loop in chiral perturbation 
theory [5] A/j = 0.35 MeV. An estimate for the upper limit of Ar including all terms of the 
order rn^ is 2 MeV [6] (it is the upper limit, because the sign of the contribution from the 
scalar meson is unknown; the conservative choice is taken in [6]). Also, to this order in the 
quark mass expansion no chiral logarithm is found. In view of this result, the approximation 
E ~ a{2y?) is good enough for the present level of accuracy. In addition to fixing the value 
of E from pion-nucleon information, one still has to consider the t-dependence of the scalar 
form factor to get an estimate for a (= a{t = 0)), so the quantity of interest is 

A,, = a{2^l') - a(0). (7) 

The first calculation was published in Ref. [7] and later more refined results were given in 
[1,5]. CHPT to one loop gives [5] A^- ~ 5 MeV, which, however, corresponds to the treatment 
of the pion-pion scattering at tree level. In [8] Acr was calculated with a dispersion relation 
with TTTT input consistent with CHPT. The result is 

A(, = 15.2 ± 0.4 MeV, (8) 

which is also obtained by Bernard et al. [9] in a calculation which includes some of the 
contributions of order q^. 

Using Eqs. 5 and 7 the a can be calculated once E is fixed from pion-nucleon information, 
namely 

(7 = E-(A^ + Ai^), (9) 

i.e. the difference between E and a is about 15 MeV with relatively small uncertainty in 
comparison with the current uncertainties in E and the CHPT result for (t(1 — y). 



DETERMINATION OF E 

Some time ago Koch used [10] hyperbolic dispersion relations to determine E and got the 
result 

E = 64±8 MeV, (10) 

where the uncertainty reflects the internal consistency of the method. To be able to estimate 
the uncertainty due to the experimental errors of the low-energy data, a dispersion method 
was proposed in Ref. [11]. That approach involves six forward dispersion relations, the stan- 
dard ones for and and, in addition, two dispersion relations for the newly defined 
amplitudes 

E^{u) = ^{A^+u:B^)\t=o, (11) 
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where a; is the pion total energy. The point here is that in the low-energy range the 2 s-waves 
and 4 p-waves matter, and the d- and f -waves can be treated as corrections (to be taken 
from an existing pliasc shift analysis; KA85). The six dispersion relations wih give six energy 
dependent functions, the partial wave amplitudes {I = 0,1). The method is a low-energy 
approach; it makes sense only to the extent that d- and higher partial waves can be treated 
as corrections. In Ref. [12] it was assumed that the d-waves were accurate to 30 %, and the 
error estimate there involved variations of the d-waves by that amount around the KA85 
values [13]. In addition to the low-energy / > 2 waves, input is needed for the higher energies. 
(The "high" energy piece starts here at the pion momentum ko = 185 MeV/c corresponding 
to a laboratory kinetic energy of 92 MeV.) In this range the imaginary parts of the invariant 
amplitudes D^, and are needed, and they are constructed from partial wave solutions 
of the Karlsruhe group (KH80, KA84 and KA85) and, also, from solutions FA93, SM95 [14] 
and SM97 of the VPI group. For the asymptotic behaviour of the amplitudes Karlsruhe forms 
[15] have been adopted. 

The six dispersion relations fully fix the amplitudes up to two subtraction constants both 
of which have their origin in the amplitude. The subtraction constant in the dispersion 
relation for £)+ is proportional to the s-wave scattering length a^_,_, and the subtraction 

constant for the E'^ amplitude to . Once these two threshold parameters are known, the 
amplitude is fully fixed. Experimental data at low energies (below the momentum ko) will be 
needed to determine a^_^_ and a^^. At the lowest energies there are results from differential 
cross section measurements which are used as input in the search. For each data set the 
contribution to the total is 



where <J*{9) is the experimental differential cross section with error /S.a{0). The computed 
cross section is denoted by cr{6), and it contains the full observable cross section constructed 
from the solutions for the hadronic amplitudes from the dispersion relations, and the elec- 
tromagnetic corrections according to the formalism of Tromborg et al. [16]. Also, the P33 
splitting between the ■n~p and ir^p channels has been taken into account. For each data set 
there is one parameter z which takes care of the normalization of that set, i.e. in the present 
analysis there is a contribution to the from the normalization. However, it does not depend 
on the number of data points in the data set. The experimental normalization uncertainty 
has been used for Az, if given in the experimental paper. The iterative procedure for solving 
the coupled dispersion relations has then the following steps: 

1. Take initial values for (aj^., a^_|_), and set c^/ = for Z < 1 and feiab < ^o- 

2. Construction of the imaginary parts of the six invariant amplitudes for all momenta. 

3. Solution of the dispersion relations. 

4. Partial waves from the real parts of the invariant amplitudes at low energy. 

5. Back to step 2 until the phase shifts stabilize. 

Usually only a few iterations are needed to get smooth partial wave amplitudes. 
It is practical to write the subthreshold expansion at the origin of the (i^, t) -plane 




(12) 



where two of the constants can be directly related to the dispersion relations 

4o = ^+(0), 4i = ^+(0). 



(13) 



(14) 
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The expression for E then gets the form 



S = 



(15) 



where the curvature term is dominated by the tttt cut giving [8] 



Ad = 11.9 ±0.6 MeV. 



(16) 



However, there is also the left-hand cut contribution, so a new ttN analysis could modify this 
number slightly. The term linear in i is a relatively sensitive quantity as can be seen from the 
numbers for the solutions A and B in [12] 



where the first figure corresponds to the dgo contribution and the second the 2^0"^ piece. 



The problem with the sigma values used to be that the value for y, the strange quark 
content, tended to be large in view of what could reasonably be expected for the strange 
quark contribution for the proton mass. Then it turned out [8,12], however, that the t- 
dependence of the scalar form factor was underestimated in the previous analyses, and that 
the new value Aa- ~ 15 MeV led to a more reasonable result, about 130 MeV of the proton 
mass coming from the strange quark piece [12]. 

This analysis relied heavily on the Karlsruhe amplitudes, which included data from the 
70's and earlier. So only a very limited amount of information from the meson factories was 
incorporated. Especially at low energies new measurements have been performed for different 
observables. New data point to some changes in the Karlsruhe amplitudes, but there is, at 
present, no general agreement of the new amplitudes, because even the new data contain 
conflicts. Also, the discussion of the value of the ttN coupling constant has continued with 
increasing vigour. In [12] it was not possible to check the dependence of the values of S on the 
pion-nucleon coupling, because one was bound to use the coupling strength corresponding to 
the input amplitudes, i.e. the Karlsruhe value /^tt = 14.28; = 0.079. The VPI group has, 
however, started to incorporate fixed-t constraints to the phase shift analysis, which gives a 
possibility to make some checks of the influence of the value of the coupling. Their value for 
the ttN coupling is = 0.076 [14]. Exploratory searches were reported in [17,18] for the VPI 
amplitudes FA93 and SM95. Definite values for S cannot be given, because the curvature 
term Ad has not been estimated with these amplitudes, but the linear part S^i gets values 
of about 3 MeV larger than the ones resulting from the use of the Karlsruhe amplitudes 
when the same low-energy data sets have been used as input. Another matter is that the 
favoured data sets would be different from the ones which are consistent with the Karlsruhe 
amplitudes. In any case, the effect of the ttN coupling seems to be relatively unimportant 
in comparison with the error bar due to other uncertainties. The reason could be that the 
S involves the l)"'"-amplitude where the Born term has been subtracted and, therefore, the 
large contribution directly proportional to is removed from the amplitude. In [18] some 
issues in the selection of the low-energy data was investigated for the case of VPI-SM95 [14] 
input amplitudes. The trend seems to be towards slightly larger values of the S^^, but the 
effect is typically about half of the estimated total uncertainty, so the overall picture does 
not change essentially. 

New low-energy data has also appeared [19,20]. The differential cross section data of Rcf. 
[19] are slightly too high in energy to be suitable for the present analysis. Rcf. [20] gives the 
flrst analyzing power results at low energy, and with more data to come [21] an important 



Ed = (-91.3 + 138.8) MeV 
Erf = (-94.5 + 144.2) MeV 



(solution A) 
(solution B) 



(17) 



RESULTS AND DISCUSSION 
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constraint on the low-energy analysis will be obtained. A careful discussion of the 7r+p data 
base has recently appeared in [22]. 

Another piece of new information has appeared quite recently, namely, a measurement of 
the backward elastic tt~p cross section [23] as a function of energy. With the KH80 input 
a good fit can be found, and the result is 'Sd = (—86.0 + 145.8) McV. This result is quite 
sensitive to the amplitudes and the P33 splitting of n^p and ir^p. Even the d-waves start to 
matter, because the VPI-SM97 d-waves differ upto a factor of two from the KA85 results. No 
attempt was made to try to estimate the error bar due to the errors in data. A nice feature 
is, however, that the resulting scattering lengths are agreeing quite well with the results from 
pionic hydrogen experiment. The search gives 

a^-p = 0.0882 fi~^ 

ao+ = 0.0922 /i'^ (18) 

to be compared with the experimental results from pionic hydrogen [24,25] 

a„-p = 0.0883 ± 0.0008 /^^^ 

0.0920 ± 0.0042 fi~^. (19) 



a, 



0+ 

Another approach for the sigma-term discussion is gradually becoming more and more 
useful. Namely, already now lattice calculations can make statements about the value of 
sigma [26,27]. Using the Feynman-Hellmann theorem 

^ dm , , 

. = »^ (20) 

the sigma can be calculated from the quark mass dependence of the proton mass. Of course, 
there are the difficulties related to the question of dynamical quarks which until now cannot 
be treated. Also, the quark mass values in actual calculations are relatively large, but gradual 
improvement can be expected here. The current numbers for a are 40-60 MeV [26], and 50 
MeV [27] where the latter calculation cites a very small error (3 MeV). 



CONCLUSIONS 

In this workshop views were expressed that the modern ttN data base is essentially inter- 
nally consistent. In the strongly constrained low-energy analysis discussed here the picture is 
not quite that clear. Also, there seems to be data sets which simply have to be "thrown away" 
in all analyses. In my opinion, one should aim at understanding from the experimental point 
of view why there is so much trouble with the n^p data. A good example is the disagreement 
of the differential cross section results with the integrated cross section measurements [28,29]. 

Another set of problems is opened by the possibility of isospin violation which has not 
been addressed here. The sensitivity of the extrapolation to the Cheng-Dashen point may 
call for an improvement of the precision in this respect. 
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